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Identifiability Conditions for Superior Detection 

S. Majidian and F. Haddadi 



Abstract — In recent years, various methods are presented for 
the case of more sources than sensors in the array processing 
which is called superior detection. Khatri-Rao approach has been 
proposed for localizing quasi stationary sources. Nested arrays 
are also used for this purpose. In this Letter, by exploiting the 
rank of covariance matrix, identifiability of these approaches 
have been discussed. Conditions for every batch and almost every 
batch of received signal has been presented. 

Index Terms — array processing, direction of arrival, identifia- 
bility, superior detection. 



I. Introduction 

DIRECTION Of Arrival (DOA) estimation for multiple 
sources has been an active research area in the past 
decades [1]. There is a vast use of Uniform Linear Array 
(ULA) in DOA Estimation. ULA is an array which elements 
are spaced equally along a straight line. Bresler in [6], shows 
limitation on the number of signals which ULA can resolve. 
This fact becomes a motivation towards the development 
of diverse classes of arrays which are Nonuniform Linear 
Arrays (NLA). Primal endeavor was by Moffet who proposed 
Minimum Redundancy Linear Arrays (MRLA) [2], which 
improves the performance of DOA estimation significantly. 
Constructing a MRLA requires an exhaustive search through 
possible combinations. In order to resolve this issue, Pal 
and Vaidyanathan presented two closed form schemes [3], 
[4], to extend the effective aperture of arrays. The so called 
nested array, which is constructed by multiple uniform linear 
subarrays (ULAs) with different inter-element spacing. The 
other scheme is formed by two sparse ULAs which satisfies 
the co-prime relationship for the interelement spacing as well 
as the number of sensors. Also Ma [5], proposed a technique 
called Khatri-Rao subspace approach, which is used for quasi- 
stationary signals. For every technique there is a crucial 
problem: the identifiability, which is the problem of possibility 
of unique direction estimation of multiple sources by an array 
of antennas. This problem is the key issues of many practical 
considerations. 

In 1986 Bresler and Mackovski [6], considered a ULA 
intercepting narrowband deterministic signals. They showed 
that the uniqueness problem is equivalent to a linear prediction 
problem and obtained a condition for a unique solution for any 
given set of observations and a condition for almost every set 
of observations. Wax and Ziskind [7], extended the results in 
[6], for any array configuration. They presented a necessary 
condition for identifiability of narrowband sources by passive 
sensor arrays. In 1992, Wax and Ziskind has studied the case 
wherein the signals are constrained [8]. 
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The problem for geolocation system is studied in [9], 
wherein multiple arrays are exploited to estimate the 2'D 
coordinates of the sources. In our model, sources are assumed 
to be in far field and therefore, their DOA is the only 
parameter of interest. Another difference is that, in [9], the 
Gaussian signal is assumed while in the Khatri-Rao approach 
investigated in this letter, this is not the case. Abramovich 
et. al. in [10], investigate the identifiability problem from a 
detection and estimation point of view. However, our analysis 
is focused on the estimation part. Because performance of the 
available methods for source enumeration such as MDL [12], 
and AIC [13], are investigated independently in the literature 
[11]. 

In this letter, we derive conditions for unique localization 
of multiple sources for Khatri-Rao, and Nested approaches. 
First, data model of signal is presented. Then, identifiability for 
every batch and almost every batch of received signal of both 
approaches are discussed. Furthermore, conditions for unique 
solution of Nested approach when there is a constraint on 
signals are presented. 

II. Khatri-Rao Subspace Approach 

A. Data Model 

Assume q narrowband signals impinging on a p-element 
ULA. Signal DOA set is 9 = {0., i = 1, ... , q}. The received 
signal x (t) is a p x 1 vector modeled as: 



x(t) = A{0)s{t) + n{t) 



(1) 



s(t), A{6) are uncorrected sources and array response 
(also known as steering matrix), respectively. n(i) is noise 
component which is assumed to be white and uncorrelated 
with source signal. In this section, we suppose that source is 
quasi stationary. It means that source variance is constant in 
each frame: 

£[|^0)| 2 ] =4 m) te [(m-l)K,mK -1] m = l,...,M 

(2) 

Where Sj(t) is i-th source. K and M are the frame length 
and the number of frames, respectively. Covariance of received 
signal x (t) in frame m shall be defined as 

R { ™ ] = E[x{t)x{t) H ] 6 C pxp Vi £ [(to - l)K, mK - 1] 

(3) 

Using signal model (1), it becomes: 



Ri m) = AR { ™ ] A H + R n 



(4) 



In which R[ m ' is the source covariance matrix in form of 

^ m) ) (5) 



Jl< m > =dmg(r^\ri m) 



<2 >• 
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A local estimation for R 



(m) 



R^. 



1 

K 



mK-l 

£ 

t=(m-l)K 



x(t)x(t) H 



(6) 



For DOA estimation in this situation, a well-known ap- 
proach called Khatri-Rao Subspace is used [5]. In this method, 
for each frame, signal covariance of received signal is vector- 
ized. 

y m 4 vec{R^) = vec(AR^A H ) + vec(R n ) (7) 
Since Ry 1 ' is diagonal, Khatri-Rao product can be used. 

y m = {A* <S A)r m + vec(R n ) (8) 
Where (•)* and ® stand for complex conjugate and Khatri-Rao 



product, respectively. r r 



_(m) (m) 



' q 



Y . Then, Y 



is constructed by collecting y m of each frame as its columns. 

Y = [ Vl ,...,y M ] (9) 
Y = (A* <g> A)ip + vec(R n )ll { (10) 



t/> = [ri,.. . ,r M ] 



r\i ri 2 
r2i r 22 



»"1M 



(11) 



Where 1^ is a row matrix of M ones. Then, noise component 
will be omitted by a projection: 

1 



P ± = Im 



AT 



MJ-M 



(12) 



At last, using Singular Value Decomposition (SVD) noise 
subspace can be used in subspace algorithms like MUSIC [1]. 



YP J 



[U s U n ] 





0" 




[Vfl 


0 


0 







(13) 



5. Identifiability For Khatri-Rao Approach 

1) Every Batch: Now condition is explained under which 
the solution will be identifiable for every batch Y. As projec- 
tion has no effect on our analyses, it is ignored. B = A* ®A. 
For supposed received signals, the DOA will be unique when: 

Y = B{9)ip T / B{e')ijj T (14) 

t.T~ 



me) b(8')} 



IT 



(15) 



If (15) holds, the nullity of \B(8) BOS')] which is assumed 
to be Q must be smaller than V = Rank 



Since 



V > i) = Rank['0 J 



(16) 



It's enough that C < rj. 

Two mutually exclusive cases for d are considered, depend- 
ing on the number of repetitions of equal directions in two 
sets of solutions: d = #pair{#j = 9j}. 

Case 1: d < 2q - (2p - 1). Since f = 2q - 
Rank[B(0) B{0')] and Rank[B(0)j is equal to minimum of 



its dimensions {q, 2p— 1}. Hence, £ = 2q— min{2p— 1, 2q—d} 
Therefore, for £ < rj, we must have: 

7] + 2p - 1 
2 



5 < 



(17) 



Case 2: d > 2q — (2p — 1). First we omit those redundant 



column where 9i = 6j and call it B 

C=(2q-d)- Rank[£?(0) B(9')] = 0 



(18) 



Since V > 0, therefore, ( < V is always true. 

2) Almost Every Batch: Now, we discuss almost every 
batch Y with exception of a set of batches of measure 
zero. Let the set of all legitimate matrices Y denoted by 
Y(6) = {Y\Y = B(0)4> T }). We assume that 0 and 9' are 
two ambiguous solution sets, we define D(9) = [j g , D(0, 0') 
where D(d,d') = {Y\Y = B{9)ij> T = B{9')ip' T }). Since 
Rank(K) = \i and every column of it is in the span of {B(9)} 
of rank 2p — 1. Therefore, (2p — complex parameters are 
required to describe Y. 



dim Y{9) = 2{2p - 1)77 
dim D(0, 9') = 2r/C 



(19) 
(20) 



Now we assume two cases for d. Case 1: d < 2q— (2p— 1). 
As = 2q — (2p — 1), dimension of D(9, 9') will be equal to 
2(2q-(2p-l))rj. Since dim (9') = q, dimension of D{0) will 
be less than or equal to 2(2q— (2p— l))r] + q. For dim D(9) < 
dim Y(9) we must have: 



q < 



4?7 



4r; + 1 



(2p - 1) 



(21) 



Case 2: d > 2q—{2p—l). In this case dimension of D(9, 9') 
will be 2drj. Therefore, dimension of D(9) will be less than 
or equal to 2dr\ + q. For dim D(9) < dim Y(0) we must 
have: 

q< 2 ^ 1 (2 P -l) (22) 

Since r/ > 1, the last inequality is more tight. 

3) Constrained Signals: In this part, we consider real 
constraint on signals. 



fk(E\s k (ti)\ : 



0 k=l,. 



,1 



(23) 



Where -B|sfc(^)| 2 is second order statistics of fcth source 
at time tj, is a smooth map from the complex plane to 
real line. A^. 1 is a vector of parameters with \i real unknowns. 
In fact, this constraint describes the relation of second order 
statistics of the sources by a function. 

The number of free parameters to describe B{0) is 2p — 
1. The source equivalent xp T is q x M real matrix with 77- 
column rank, which is described by qrj free real parameters. 
Let the set of all tp, rank-?? matrix which comply with above 
constraint, be denoted by D. The number of free parameters is 
reduced by constraint to q-q — qp,. Therefore, dimension of D 
is qrj — {qrj — qp) = q\i. Set of ambiguous solution is similar 
to "almost every batch" section defined DOS) = [j 0l . 

For identifiability, its enough that the dimension of the set 
D{6) be smaller than dimension of D. Similar to the last 
section, in two cases we have: 
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Case 1: d < 2q - (2p - 1). As ( = nuU[£(0) B(8% 
and Raak.[B(8)B(6')] = min{2p - l,2q - d}. Hence, 
Q = 2q — (2p — 1). Thus, 2(2q — (2p — 1)) real parameters 
are required to describe an arbitrary vector in null space of 
[Bid) B(0')] and 2(2g- (2p-l))r? are required to describe 



an arbitrary rank-77 matrix 



/'/" 



in this null space. Therefore, 



dim D(9, 9') = 2(2q—(2p—l))rj—(q-q — qp). Since dimension 
of D{6') is q, Therefore 

dim D(8, 9') = 2{2q - (2p - l))r? - (qrj - qp) (24) 

For dim D(8) < dim D we must have : 

2r? 



9 < 



3r; + l 



(2p - 1) 



(25) 



As we can see, number of unknown parameters of constraint 
have the same effect on dimension of d and dim D(9). While 
it does not affect the condition. 

Case 2: d > 2q — (2p — 1). According to the last section, 
we omit the redundant columns. Therefore 



(=(2q-d)- Rank[B(0) B(8')\ = 0 



(26) 



and dim D(8, 9') = 0 - (qrj - qp). Therefore, dim D{8) < 
-{qr\ - qfi) + q. If r] > 1 then dim D{8) < dim D. 

4) Alternative Proof For Constrained Signals: ip T is de- 
scribed by qr] free real parameters. The number of free param- 
eters is reduced by constraint and becomes (qrj — qp). Since 
xj) is affected by q source, the dimensionality of legitimate set 
is equal to qp + q ambiguous solution: 



C = B(9)ip T ^ B(9')xp 1 



The ambiguity set is defined by 



D = {B{9)x)) T : C = 0 constrained} 



(27) 



(28) 



B(9)ip T has 2qr) + q complex parameters. The constraint 
reduces it by (2p — 1)77. Therefore, its dimension will be 
2(2qr]—(2p— l)rj + q). For unique solution, the dimensionality 
of the ambiguity set should be smaller than dimensionality of 
legitimate set, which means 2{2qrj — (2p — 1)77 + q) < qr\ + q 
which results in (25). 



III. NESTED APPROACH 

A. Data Model 

Two-level nested array consists of two ULA, in which the 
inner one has L\ elements with spacing d\, and the outer has 
L2 with spacing di = (Li + l)di. Therefore, the location of 
each element is 



d = [di,2di, . . ..Lxdx, (Li + l)di, 

2(Li + l)d 1 ,...,L 2 (L 1 +1)^] 



Define 



D = {di - dj\i,j = 1, . .. ,L = L\ + L 2 } 



This array can act as a virtual array with this element: 

d v = [-kdi,..., kdi] k = L 2 (L 1 + l)-l 



(29) 
(30) 
(31) 



In the nested approach, covariance matrix of received signal 
is 

R x = AR S A H + a 2 I L (32) 

Where a 2 is the noise variance and Rs = diag(erf , . . . , a 2 ) 
is covariance of uncorrected sources. Now, z is defined as: 



Consequently 



z = \ec(R x ) 



z = (A* ®A H )p + cj 2 l L 



In which p = [a 2 ,..., a 2 ] , and 1^ 



(33) 



(34) 



iff 



where e, is a vector with the i-th element of 1 and the others 
of 0. Construct A\ by omitting redundant rows of A* 
Therefore 

zi = Aip + a 2 e' 



(35) 



Where e' is a (L 2 + 2L-2)/2 x q vector in which the L 2 /4 + 
L/2-th element is 1 and the others are 0. 

B. Identifiability For Nested Approach 

1 ) Every Batch: The condition will be derived under which 
the solution will be unique for every batch of data. 

z 1 = A 1 {8)p^A' l {8') V f 



[AM A[(8')} 



7^0 



(36) 
(37) 



If (37) holds, nullity of [Ai(8) A'^8')] (= 0 must be 



smaller than V = Rank 



Since V > rj = Rank[p], It is 

enough that Q < 77 also in this situation, we have two mutual 
exclusive cases for d. For case d < 2q — (L 2 + 2L — 2)/2: 



Q < 



rj+{L 2 + 2L-2)/2 



(38) 



For the complement case, it can be shown that, there is no 
restrictive condition. 

2) Almost Every Batch: It is straightforward to show: 

2r) 



'/< 2 ?? Tl (L2 + 2L_2)/2) 



3) Constrained Signals: 
2'i] 



9< 3^TT (i +2L " 2)/2) 



(39) 



(40) 



IV. Conclusion 

In this letter, conditions for identifiability of signals in both 
Khatri-Rao and Nested approaches for every batch, almost 
every batch, and the constrained data were derived. It is proved 
that the number of resolvable sources is 0(L 2 ) and 0(2p) 
for Nested and Khatri-Rao approaches, respectively. The gap 
between the performance of the existing methods and the 
upper bounds in this letter can be closed by development 
of better algorithms for superior detection in array signal 
processing. 
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